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Abstract 

Fully differential heavy quark contributions to the photon structure functions 
in deeply inelastic scattering are computed in next-to-leading order QCD, in- 
cluding both the direct and resolved contributions. A variety of distributions 
are presented and discussed. Several of the distributions show marked differ- 
ences between the resolved and direct cases due primarily to the presense of 
the gluon distribution in the former and lack thereof in the later. 
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I. INTRODUCTION 



Hard scattering processes in two photon collisions have long been recognized as playing an 
important role in the investigation of QCD. Examples include the photon structure function, 
the production of high-p^- jets, and heavy quark production, to name just a few. Such 
processes can yield information pertaining to the hard-scattering parton-level cross sections, 
the distributions of partons in one or both initial photons, and the fragmentation of partons 
into hadrons in the final state. The information obtained thereby is complementary to that 
obtained in other hard scattering processes. New developments in these and related areas 
have recently been reviewed in ||T|. 

Activity in the area of two photon physics has increased recently with new data becoming 
available from several groups and with the initiation of the LEP2 research program. In 
particular, it is anticipated that substantial samples of charm quark production data will 
soon be obtained. This data sample can be divided into three subsets, each of which is 
treated differently from a theoretical standpoint. First, both photons which originate from 
the e"*" and e~ beams can be nearly on-shell. This is by far the largest subset and the relevant 
theoretical issues have been discussed in 0, for example. The second largest subset, and 
the one addressed in this paper, consists of one photon being nearly real and the other being 
substantially off-shell, with a negative squared invariant mass satisfying > 2 (GeV/c)^, 
for example. The third subset consists of both photons being far off-shell, but the rate is 
substantially reduced and this case will not be addressed here. 

When an on-shell photon interacts in a hard-scattering process, two cases must be con- 
sidered. First, the photon may interact via the pointlike electromagnetic coupling to quarks. 
This component is often referred to in the literature as the direct or pointlike piece. In the 
case considered here, this contribution is calculable in perturbation theory once the heavy 
quark mass and the strong coupling are specified. Second, the photon may undergo a fluc- 
tuation into a hadronic state. The connection to the hard scattering is then made through 
parton distributions in the photon. These can be obtained, for example, from studies of 
the photon structure function or the production of high-p^ jets in two photon collisions or 
in photoproduction. On the other hand, the sensitivity of heavy quark production to these 
distributions may be used to provide additional constraints on them. 

The purpose of this work is to present the fully differential heavy quark contributions 
to the photon structure functions. These may be compared directly with experimentally 
measured structure functions, should they become available, or embedded in a larger Monte 
Carlo program for generation of heavy hadron events. Heavy quark correlations have also 
been studied theoretically in hadroproduction P], photoproduction 0, electroproduction 
1^, and recently in direct photon-photon collisions |^ as well as experimentally 

Our calculation is performed using the subtraction method which is based on the replace- 
ment of divergent (soft or coUinear) terms in the squared matrix elements by generalized 
distributions. The method was first used in the context of electron-positron annihilation |^ 
and its essence is nicely described and compared to the phase-space slicing method in 
the introduction of a paper by Kunszt and Soper [|1^ . 

In Sec. II the method of calculation is reviewed and analytical results presented. Nu- 
merical results and a discussion of related physics issues are presented in Sec. HI and some 
conclusions are given in Sec. IV. 
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II. METHOD 



Deep inelastic photon structure functions F^lxjQ"^) {k = 2,L) are measured in e^e 
collisions via the process 

e-(/)+e+^e-(/') + e+ + X (2.1) 

where X denotes any hadronic final state allowed by quantum number conservation. See 
Figure 1. When the outgoing electron is tagged the reaction is dominated by 

Y{q)+l{p)^X (2.2) 

where one of the photons is highly virtual and the other is near mass shell. For heavy quark 
production the relevant reaction is 

7*(g) +7(P) Q{Pi) + QiP2)+X (2.3) 
and the cross section is given by 
dV(e+e^ QQe+e~X) ^-na^ 



dxdy 



-S 



{l + (1 - yf} F^{x, Q\ m') - y'F2{x, Q\ m^)] (2.4) 



after integration over the azimuthal angle between the plane containing the incoming and 
outgoing electron and the plane containing 7(p) and the outgoing heavy quark. Here q = 
I — I', = —q"^, X = Q'^/2p ■ q, and y = p ■ q/p ■ I- The heavy quark contributions to the 
photon structure functions are denoted by F^{x,Q'^,m'^) {k = 2,L) with the heavy quark 
mass m shown explicitly. The electromagnetic coupling is a = e^/47r and S is the square of 
the c. m. energy of the electron-positron system. 



We follow the terminology of [|ri|] in distinguishing between direct and resolved photon 
contributions. Direct contributions are those reactions where the incident photon partic- 
ipates directly in the hard scattering process. Resolved contributions are those reactions 
where the incident photon participates in the hard scattering via one of its constituents. 
For heavy quark production at next-to-leading order, the direct component is independent 
of the parton distributions in the photon. This is expressed as 

Filix, Q\ m\ /ij) = Fl^'ix, Q\ m\ ^)) + F,^'^(x, Q^ m^). (2.5) 

The mass factorization scale /xj is explicitly shown. It is absent in the second term because 
any potential coUinear singularities are regulated by the heavy quark mass m. Calculation 
of Fl'^ and F^'^ in fully differential form at NLO will be discussed in the following sections 
and represents a generalization of the fully and single inclusive results already available in 

In our computation we use the subtraction method which is based on the replacement 
of divergent (soft or coUinear) terms in the squared matrix elements by generalized dis- 
tributions. This allows us to isolate the soft and coUinear poles within the framework of 
dimensional regularization without calculating all the phase space integrals in a space-time 
dimension n 7^ 4 as is usually required in a traditional single particle inclusive computation. 
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In this method the expressions for the squared matrix elements in the colhnear hmit appear 
in a factorized form, where poles in n— 4 multiply splitting functions and lower order squared 
matrix elements. The cancellation of colhnear singularities is then performed using mass 



factorization [|13|. The expressions for the squared matrix elements in the soft limit also ap- 
pear in a factorized form where poles inn — A multiply lower order squared matrix elements. 
The cancellation of soft singularities takes place upon adding the contributions from the 
renormalized virtual diagrams. Since the final result is in four-dimensional space time, we 
can compute all relevant phase space integrations using standard Monte Carlo integration 
techniques and produce histograms for exclusive, semi-inclusive, or inclusive quantities 
related to any of the outgoing particles. We can also incorporate various experimental cuts 
defined in terms of partonic variables. 



A. Direct component 

To order as, the subprocesses where the incident on-shell photon participates directly in 
the hard scattering process are given by 

7*(g) 7(^1) ^ Q{pi) Q{P2) (2.6) 

7*(g) 7(fci) ^ Q{pi) Q{P2) g{h). (2.7) 

For future reference we introduce the invariants s = (g + kiY and S5 = {pi +P2Y ^'^d define 
additionally, s' = s — and = — q^. The direct piece of the fully differential heavy 
quark contribution to the photon structure function follows from 

dF2'''{x,Q\m^) = ^da, (2.8) 

with k = 2,L and da2 = dac + SdaL/^. The virtual photon-photon cross sections are 
summarized as 

dai = daf^ + da^"^ + da\^^ (2.9) 



for i = G, L. The first term is the lowest order contribution from (|2.6| ) and is given by 

daP = Ci Mf^^ dT2 (2.10) 

where Cq = l/{^s') and Cl = l/(4s'). dV2 is the two body differential phase space factor. 
m/^-* = 1287r^ ej^a'^NBi QED where en is the charge of the heavy quark in units of e, = 3 
for SU{3), and Bi qed may be found in 0, for example. 

The second term in ( p.9|) is the finite contribution remaining after the addition of the 
soft contribution from ( p.7| ) and the interference of ( |2.6| ) with its renormalized one loop 
corrections. It is expressible as 

rfaf") = Ci Mt^ dT2 (2.11) 
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where Mj*-™'' = Airascj^a'^NCFVi^QED- cts is the strong couphng and Cp = 4/3 for SU{3). 
The expression for Vi^QEo is too long to be presented here but may be obtained from the 
authors upon request. The third term is a finite contribution given by 

da^f^ = 2(^)'Qf3,j^(^) Mh^lfll^dzdy sin e.de^de^. (2.12) 
V167r2y *'"5(s')2 Vl-2/p(l + 2/)(l-y) y 1 1 2 K J 



Here we have defined (3^ = \Jl — Am? / s^^ z = sy s', and y denotes the cosine of the angle 
between q and k2 in the 7*7 center of mass frame. fi{z, y, 9i, 62) is the coefficient of C% in 
/f (2;, 6'i, 6*2) given in it is, up to overall factors, the squared matrix element for the 
subprocess in Eq. (|2.7|) times (1 — zY{\ — y)(l + y). Note that fi is finite in the soft {k,2 0) 



limit. The p distribution is defined by 



f ^) ./'../M^. (,13, 



p \L — p J p 



The lower limit of the RHS of (p.l3|) , when inserted in (|2.12|) , cancels exactly against a 



corresponding piece in Vi^QED so that the overall expression is independent of p. 



B. Resolved component 



To order a^, the subprocesses where the incident on-shell photon participates in the hard 
scattering via one of its constituents are 



7 9 

* 

7 9 

* 

7 9 

* — 

7 (I 



QQ 
QQg 
QQq 
QQq 



(2.14) 



The calculation of the fully differential heavy quark contributions to the photon structure 
function proceeds along the same lines as those for the proton structure functions and will 
not be repeated here. The interested reader is invited to consult for full details. The 
final result is of the form 



^,^'^(x,Q^m^/i2) 



Q^asip 





i=q,q 



k,i 



+ efd^] 



I (1) 



(2.15) 



where k = 2,L. The lower boundary on the integration is ^mm = x{Am? + Q'^)/Q'^- The 
parton momentum distributions in the photon are denoted by fi/-y{^, p"^)- The mass fac- 
torization scale pf has been set equal to the renormalization scale pr and is denoted by p. 
Finally, cj^\ and cj^j , (/ = 0, 1), and d^f}] and o^\- are scale independent parton coefficient 
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functions. In Eq. (|2.15|) the coefficient functions are distinguished by their origin: Cj^\ and 

'cj!'l originate from the virtual photon-heavy quark couphng and therefore appear for both 

charged and neutral parton-induced reactions; arise from the virtual photon-light quark 

coupling; o[.\- come from the interference between these processes. All charges are in units 

of e. We have isolated the scale dependent terms, proportional to c[.\^ Finally, note that 
Eq. ( |2.15| ) only holds for Q"^ > 0. When = one encounters double resolved processes 
which are not discussed here. 



III. RESULTS 

Using the results of the previous section we have constructed a program for the fully 
differential heavy quark contributions to the photon structure functions. With this program 
we can calculate fully inclusive, single inclusive, and fully differential structure functions. 
The program uses Monte Carlo integration so it is possible to implement experimental cuts, 
provided that they are defined in terms of partonic variables. Additional integrations and 
final state fragmentation may be added to produce hadronic cross sections. We have verified, 
to better that 1%, that our program produces fully inclusive structure functions that match 



those already available |T5 



We use a mass of 1.5 GeV for the charm contribution to the photon structure functions 
and neglect the bottom quark content. The latter contribution is suppressed by the larger 
mass and smaller charge of the 6— quark. The strong coupling as evaluated with Nf = 3 
at one loop accuracy is used in the LO results. For the NLO results the two loop value is 
used. The GRV photon parton distribution functions |]TB|] in the MS scheme are used. This 
parameterization has A^^"^'' = 232 MeV and A^"^"* = 248 MeV for the leading order and 
higher order parameterizations, respectively. 

In Fig. 2 the direct and resolved contributions to F2{x,Q'^,ml)/a at = 10 (GeV/c)^ 
are shown as functions of x. We choose the renormalization and factorization scales to be 
equal to fiQ where fi^ = + 'imi. This choice reduces to yU^ = Q"^ for the electroproduction 
of massless quarks and /i^ = 4m^ for the photoproduction of charm quarks. The dashed 
lines are the LO contributions and the solid lines are the NLO contributions. Clearly, the 



two components dominate in different x regions, as was noted previously in [|I2| . The direct 
component peaks just above threshold while the resolved component peaks at small x due 
to a sharp rise in the gluon density in the photon. 

Typically, the photon-quark initiated process contribution to F2{x,Q'^,ml) is so small 
that it may be neglected. However, this is a scale dependent statement which must be veri- 
fied. Using our results we have checked that formally the variation of F2{x,Q^,ml, fi^, fi'j) 
with respect to /i^ or /i/ is one higher power in Ug- We have also studied this numerically by 
keeping fir and fif separate and distinct and examining the change in the result under their 
variations. At x = 10"^, when /x/ is varied between 2/io and /io/2, -^^(x, Q^, m^, /i^ = /Iq, /ij) 
changes by +6% and —14%, respectively, at NLO to be compared with +26% and —29% 
at LO. Under the same variation, the sum of quark and anti-quark initiated process contri- 
butions changes from —10% to +5%. These contributions are around —1% at the central 
value. Also at a; = 10"'^ when /i^ is varied between 2/io and /io/2, -F^ (x, Q^, ml, fif., fi'j = /Xq) 
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changes by —11% and +13%. respectively, at NLO to be compared with —19% and +31% 
at LO. Under the same variation, the sum of quark and anti-quark initiated process contri- 
butions changes from —1% to —2%. We therefore conclude that at small x the dominance of 
the resolved over the direct contribution to F2 (x, Q^, ml) could provide information on the 
gluon content of the photon with no need to consider the (anti)quark initiated processes. 

Let us next examine a sample of the various differential structure functions that can be 
produced with our resuhs. For these plots we have not applied any cuts. We choose the 
renormalization and factorization scales to be equal to iiq with Hq = + 4m^ + where 

= {Pt{c) + P^{c))/2. Note that there are many possible choices of scale. We have simply 
chosen to add in the average of the squares of the transverse momenta of the heavy quarks. 
For illustration purposes, we consider the limiting cases of small and large x in the following 
figures. In all plots we take the z axis to be in the direction of the 7* in the 7*7 c. m. frame. 

We begin by examining the invariant mass distribution of the charm- ant icharm pair. 
The results for d{F^{x, Q^, ml)/a)/dMcc at = 10 (GeV/c)^ are shown in Fig. 3 with the 
LO results shown as dashed lines and the NLO results shown as sohd fines. At x — 10~^ 
the resolved contribution dominates whereas at a; = 0.3 the direct contribution dominates. 
For the resolved case, the NLO curve is flatter than the LO curve because the quark and 
anti-quark initiated processes reach their largest negative value just above threshold. This 
compensates for the gluon initiated process which tends to make the curve steeper due to 
negative contributions at large invariant mass. For the direct case, at leading order, the 
invariant mass is fixed at ^/s with s — (5^(1 — x)/x. At NLO, the radiated gluon prefers to 
be soft so the correction rises from zero at threshold 2mc to a maximum at ^Q'^{1 — x)/x. 

The transverse momentum distributions (i(F2'(a;, Q^, m^)/a)/(iP((c) at = 10 (GeV/c)^ 
are shown in Fig. 4. Pt{c) is the transverse momentum of the charm quark relative to the 
z axis defined above. The LO results are shown as dashed lines and the NLO results are 
shown as sofid fines. The results for x — 10~^ are shown in Fig. 4(a) and for x — 0.3 in 
Fig. 4(b). Also shown for comparison in the inset is the LO result for Q'^ — 100 (GeV/c)^ 
and X = 0.01. For the direct case there is a Jacobian peak at Pt{c)^^^ = \J s/A — ml and a 
factor of Pt{c) which forces a zero at Pt{c) = 0. Depending on a peak may develop also 
at small but non-zero Pt{c) as shown in the inset. However, at these values oi x and the 
resolved contribution dominates. 

At leading order d{F2{x,Q'^,ml)/a)/dPt{cc), shown in Fig. 5 at = 10 (GeV/c)^, is 
a delta function at Pt{cc) — for both the direct and resolved cases. Pt{cc) is the sum of 
the charm and anticharm transverse momenta and is identically zero at LO but is balanced 
by the transverse momentum of the radiated gluon at NLO. The resolved contribution 
has large negative contributions at Pt{cc) = due to a combination of soft and collinear 
subtractions. This accounts for the reduced value in the first bin. If fact, for very narrow 
binning these distributions become negative at Pt{cc) — 0. The direct contribution has only 
soft subtractions which are much less extreme. 

Another distribution that shows similar effects is d{F2 {x , Q"^ , m"^) / a) / d{A(f)) shown in 
Fig. 6 at = 10 (GeV/c)^. A0 is the azimuthal angle between the charm and anti-charm 
in the plane perpendicular to the z axis. At LO, for both the resolved and direct cases, 
this distribution is a delta function at the back-to-back configuration A0 = tt. Again there 
is a depletion at NLO in the bin nearest A0 = tt. This is most noticeable in the resolved 
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contribution where there are both coUinear and soft subtractions. 

The role of parton distributions in the photon can be further seen by examining 
various (pseudo) rapidity distributions. d{F2{x,Q'^,m'l)/a)/dy{c) is shown in Fig. 7 and 
d{F2'{x,Q'^,ml)/a)/dr){c) is shown in Fig. 8, both at Q'^ — 10 (GeV/c)^. The charm quark 
rapidity y{c) is defined by y = l/21n[(£' + pz)/{E — pz)] and the pseudo-rapidity ri{c) is 
defined by 77(c) = l/21n[(l + cos^)/(l — cos^)]. AU quantities arc measured relative to the 
z axis in the 7*7 c. m. frame. The dominance of the gluon distribution at small x together 
with the LO matrix element for ■j*g — > QQ produces a resolved distribution for y{c) shown 
in Fig. 7(a) that is offset from zero. In contrast, the absence of the gluon distribution in 
the direct case, shown in Fig. 7(b), gives a relatively flat distribution. The comparison of 
these two distributions offers a striking illustration of the effect of the gluon distribution in 
the photon. Data for these distributions may well play a significant role in constraining the 
gluon distribution. 

Changes in the above differential distributions with variations in the scale are of two 
types. First, if the distribution is trivial at LO (i.e. delta function hke) then the effect of a 
scale variation is an overall shift in the normalization of the curve with little effect on the 
shape. However, if the distribution is non-trivial at LO, then the effect of a scale variation 
is mild changes in both shape and normalization. 

IV. CONCLUSION 

In this paper we have presented the fully differential heavy quark contributions to the 
photon structure functions of deeply inelastic scattering in next-to-leading order QCD. The 
calculation was performed using the subtraction method. A variety of distributions were 
discussed for both the resolved and direct components. Upon integration, our fully differ- 
ential results reduce to, and agree with, the single inclusive results already available in the 
literature. 

We studied the renormalization and factorization scale dependence of our results. At the 
fully inclusive level, the scale variation of F2 {x , Q'^ , m'^) is greatly reduced at NLO relative 
to the LO calculation. In addition, the sum of the quark and anti-quark initiated process 
contributions was shown to be small. 

The distributions shown in this paper show rather striking differences for the direct 
and resolved cases. This is due primarily to the presence of the gluon distribution in the 
resolved contribution and lack thereof in the direct contribution. It may therefore be possible 
to obtain information about the gluon distribution in the photon that is complimentary to 
that obtained in other hard scattering processes. 
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Figure Captions 
FIG.l. The process e"*" + e~ — e+ + + X. 

FIG. 2. The direct and resolved contributions to F^{x,Q'^,ml)/a at = 10 (GeV/c)^ 
as a function of x. The direct contribution dominates at large x and the resolved 
contribution dominates at small x. The dashed lines are the LO results and the solid 
hnes are the NLO results. 

FIG. 3. The contributions to d{F^{x,Q^,ml)/a)/dM^ at = 10 (GeV/c)^. a.) Predomi- 
nantly resolved at x = 10~^ and b.) predominantly direct at x = 0.3. The dashed lines 
are the LO results and the sohd hnes are the NLO results. 

FIG.4. The contributions to d{F^{x,Q^,ml)/a)/dPt{c) at = iq (GeV/c)^. a.) Predom- 
inantly resolved at a; = 10~^ and b.) predominantly direct at x — 0.3. The dashed 
hnes are the LO results and the solid hnes are the NLO results. The inset shows the 
LO direct contribution evaluated at Q"^ — 100 (GeV/c)^ and x — 0.01. 

FIG. 5. The contributions to d{F;'{x,Q'^,ml)/a)/dPt{cc) at Q'^ = 10 {GeV/cf. a.) Pre- 
dominantly resolved at x — 10~^ and b.) predominantly direct at x — 0.3. The dashed 
hnes are the LO results and the solid hnes are the NLO results. 

FIG.6. The contributions to d{F^{x, Q^, ml)/a)/d{A(p) at = 10 {GeV/cf. a.) Predom- 
inantly resolved at x — 10"^ and b.) predominantly direct at x — 0.3. The dashed 
lines are the LO results and the solid hnes are the NLO results. 

FIG. 7. The contributions to d{F^{x,Q'^,ml)/a)/dy{c) at = 10 {GeV/cf. a.) Predomi- 
nantly resolved at a; = 10~^ and b.) predominantly direct at x — 0.3. The dashed lines 
are the LO results and the solid lines are the NLO results. 

FIG.8. The contributions to d{F^{x,Q^,ml)/a)/dr]{c) at = 10 (GeV/c)^. a.) Predom- 
inantly resolved at a; = 10~^ and b.) predominantly direct at x — 0.3. The dashed 
lines are the LO results and the solid hnes are the NLO results. 
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